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Abstract 

The linear singularly perturbed convection-diffusion problem in one 
dimension is considered and its discretization on a Bakhvalov-type mesh 
is analyzed. The preconditioning technique is used to obtain the pointwise 
convergence uniform in the perturbation parameter. 

Keywords: singular perturbation, convection-diffusion, boundary-value prob¬ 
lem, Bakhvalov-type mesh, finite differences, uniform convergence, precondi¬ 
tioning 

2000 MSC: 65L10, 65L12, 65L20, 65L70 

1 Introduction 

The report is a supplement to [5]. 

2 The continuous problem 

We consider the problem 

Cu := —eu" — b{x)u' + c(x)u = f{x), x £ (0,1), u{0) = u{l) = 0, (1) 

with a small positive perturbation parameter e and (^^[0, l]-functions b, c, and 
/, where b and c satisfy 

b{x) > /3 > 0, c{x) >0 for X G / := [0,1]. 

*This author’s work is supported by the Irish Research Council under Grant No. 
RS/2011/179. 
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It is well known, see B 0 for instance, that (ED has a unique solution u in 
C^{I), which in general has a boundary layer near x = 0. Our goal is to find 
this solution numerically. 

The solution u can be decomposed into the smooth and boundary-layer parts. 
We present here Linfi’s [H Theorem 3.48] version of such a decomposition: 

u{x) = s{x)+y{x), ( 2 ) 

\s^'‘\x)\ < C {1 + , |j/W(a:)| < (3) 

X G I, k = 0,1,2, 3. 

Above and throughout the report, C denotes a generic positive constant which 
is independent of £. For the construction of the function s, see [3], since the 
details are not of interest here. As for y, it is important to note that it solves 
the problem 


Cy{x)=0, a:e(0,l), y(0) =-s(0), ?/(l) = 0, 

with a homogeneous differential equation. We shall use this fact later on in the 
report. 


3 The discrete problem and condition number 
estimate 


We first define a finite-difference discretization of the problem ED on a general 
mesh with mesh points Xi, i = 0,1, ..., N, such that 0 = xq < xi < • • • < 
xn = 1- Throughout the rest of the paper, the constants C are also independent 
of A^. 

Let hi = Xi — Xi-i, i = 1,2,... ,N, and hi = {hi-\-hij^i)/2, i = 1,2,..., N—1. 
Mesh functions on are denoted by , etc. If 5 is a function defined 

on I, we write gi instead of g{xi) and g^ for the corresponding mesh function. 
Any mesh function is identified with an (A^-l-l)-dimensional column vector, 
= \Wq , Wi ,... , and its maximum norm is given by 

IIW^II = max jW^^j. 

0<i<Af 


For the matrix norm, which we also denote by || • jj, we take the norm subordinate 
to the above maximum vector norm. 

We discretize the problem ED on using the upwind finite-difference 
scheme: 

= 0 , 

■.= -£D''U^ -hD'Uf + c,U^ = i=l,2,...,N -1, (4) 

U^ = 0, 


where 


and 




1 


’i+i 

hi 


i+1 


hi 


D'Wl^ = 

hi+i 
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The linear system (jT]) can be written down in matrix form, 


AnU^ 


= f 


-N 


( 5 ) 


where An = [aij] is a tridiagonal matrix with ooo = 1 and unn = 1 being 
the only nonzero elements in the 0 th and A'^th rows, respectively, and where 
/'^ = [0,/l,/2,...,/iV-l,0]^. 

It is easy to see that An is an L-matrix, i.e., an > 0 and aij < 0 if i j, for 
alH, j = 0,1,..., -/V. The matrix An is also inverse monotone, which means that 
it is non-singular and that A'^ > 0 (inequalities involving matrices and vectors 
should be understood component-wise), and therefore an M-matrix (inverse 
monotone L-matrix). This can be proved using the following M-criterion, see 
[ 2 ] for instance. 

Theorem 1. Let A be an L-matrix and let there exist a vector w such that 
w > 0 and Aw > 7 for some positive constant j. A is then an M-matrix and it 
holds that ||^“^|| < 7 “^||w||. 

To see that An is an M-matrix, just set Wi = 2 — Xi, i = 0,1,..., iV in 
Theorem [T] to get that Anw > min{l,/3}. This also implies that the discrete 
problem ([S]) is stable uniformly in s, 






< C. 


( 6 ) 


Of course, the system has a unique solution . 


4 A Bakhvalov-type mesh 

A generalization of the Bakhvalov mesh [T] to a class of Bakhvalov-type meshes 
can be found in [9] . Here we take one of the Bakhvalov-type meshes from [9] for 
the discretization mesh . We refer to this mesh as as Vulanovic-Bakhvalov 
mesh (VB-mesh). The points of the VB-mesh are generated by the function A 
in the sense that Xi = X{ti), where L = i/N. The mesh-generating function A 
is defined as follows: 


X{t) = 


tpia) -\- ip'{a){t 


t G [0,a], 
a), t G [a, 1], 


(7) 


with 0 < g < 1 and if = aef), where 

(j){t) = ^ - 1, t€ [0, a], 

q — t q — t 

On the interval [a, 1], A is the tangent line from the point (1,1) to ip, touching 
Ip at {a,ip{a)). The point a can be determined from the equation 


ip{a) -\- tp'{a){l — a) = 1. 

Since (p'{t) = q/{q — t)'^, the above equation reduces to a quadratic one, 
aea{q — a) -\- aeq{l — a) = {q — a)^. 
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which is easy to solve for a: 


a = 


q — \/ aeq{l — q + ae) 
1 “t" CLE 


We have to assume that ae < q (which is equivalent to < 1) then 

a > 0. Note also that a < q and 


9-a = C\/e, C<C, 




( 8 ) 


Let J be the index such that tj_i < a < tj. Starting from the mesh point xj, 
the mesh is uniform, with step size H. However, xj behaves differently from 
the transition point of the Shishkin mesh because 


xj > ipia) 



We note that the transition point is different also from the Bakhvalov- 

Shishkin of Vulanovic-Shishkin meshes in the sense of [7]. 

We now give the estimate for the condition number of An when the dis¬ 
crete problem (|3]) is formed on the VB-mesh as described above. The condition 
number is 

k{An) := ||AjY^||||HAr||. 

We estimate the upper bound for ||Ajv|| by examining the entries of the matrix 
An directly, 

Ar2 

||A^||<C—. 

£ 

Combining this with ([S]), we get the following result. 

Theorem 2. The condition number of An on the VB-mesh satisfies the follow¬ 
ing sharp bound: 

N'^ 

k(An) ^ C —. 

£ 


5 Conditioning 


Let M = diag (ruo, toi, ..., toat) be a diagonal matrix with the entries 


Too = 1, 



i = 1, 2,..., iV — 1, and rriAr = 1. 


In other words, 


mo = 1, 



j = 1, 2,..., J, and m^ = 1, z = J -|- 1,..., iV. 


(9) 


When the system m 
equations 1,2, ..., J 
modified system is 


is multiplied by M, this is equivalent to multiplying the 
of the discrete problem (|3]) by hi/H, i = 1, 2,..., J. The 

AnU^ = Mf^, ( 10 ) 
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where Aj^i = MA^. Let the entries of Ajv be denoted by Uij, the nonzero ones 
being 


H ■— j ^ 


h,H' 

l<i< J-1, 

S 

hjH^ 

i = J, 

€ 

J + l<i<N-l 


bihi 


hiA-^H hi+iH 


Ti .— — < 


H+1 

_s _ 

_e _ ^ 

m H’ 


and 


di •— da — ^ 


1 , 


i = J, 

J + l<i< N -1, 
i = 0 


hi 


h Ci^ 1 ^ ^ ^ 

—h — + J 1 ^ i ^ JV — 1, 


L 


i = N. 


Unlike the Shishkin mesh, which is piece-wise uniform, the VB-mesh is 
graded in the fine part. Because of this, it is more difficulty to prove the uniform 
stability of the modified scheme. This is done in Lemma [5] below, but first we 
need some crucial estimates for the graded mesh defined by ©• 

Lemma 1. For the mesh-generating function given in m, the following esti¬ 
mates hold true: 

e{hi+i - htj_ ^2 


hih 




and 


hjH 


f = l,2,...,J-2, 

(11) 

VI 

(12) 


Proof. For i < J — 2, we have 
hi = Xi — Xi-i = ae 

hi+i = 


q-ti q- ti-i 
aeq 


aeq 


N{q - U-i){q - UY 


and 


hi+i hi — 


N{q-ti){q-U+i)' 

2aeq 


N‘^{q - ti_i){q - ti){q - ti+i)' 
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Then (fTTl) follows because 


e{hi+i-hi) 2 {q-ti) 


hjh 




aq 




The proof of (fT^ is more complicated due to the presence of hj. First, 
hj = + 72, where 71 = Xq, — xj-i, 72 = iCj — Xa, and Xa = ' 0 (a). Since 

72 = ip'{a){tj - a) 
aeq f tj — a 


and 


a \ q — a 


71 = ae (0(a) - 0(tj_i)) 
/ a t,/_i 


= as ' 


we have 


hj = 



' \q- 

-a g - tj- 


aeq 

a — tj-i 


q — a 

q-tj-i ’ 

aeq 

tj - a 

a — tj-i 


q — a 
aeq 


q — a 
tj — a -\- 


q-tj-i 

{q - a)(a - tj i) 


(g- a)2 

tj -a + 


aeq 


q - tj-i 
cy£(a-tj_i) ' 
9-ij-i 


Moreover, 


0'(a) = 




(g - ay 


and H = xj+i — xj = 




N 


implying that 


H = 


aeq 


N{q— ay 


Therefore, 


H-hj = 


1 


tj — a a — tj-i 


q - tj-i 
1 


aeq 
q — a 
aeq 
q — a 
aeq 
q — a 

asq , . X 

(a tj—iy, ,, , 

q-a (g - a)(g - tj_i) 

aeq (a-tj_i)^ 


N{q-a) q-a q-tj-i 
a — tj_i a — tj-i 
q — a 

(a - tj-i) 


1 


q-a g-tj_i 
a — tj-\ 


(g - a)2 


q - tj-i 
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We now have 


H — hj ae^q {a — tj-if' q — a 1 (q — a)'^N 

hjH {q — q — tj-i aeq t^-a _|_ a-tj-i 

Q Ot- Q tj—l 

_ {q - a)N ^ {a - _ (g - a){q - tj-i) 

aq q-tj^i + 2atj-i - 

_ {q-afN ^ (g - tj-if ^ (^e ^ 

aq oj ~ aqN uj ’ 


where 


w := ^ - g^ + 2 gtj_i - tj-itj 

and where in the last step we used ([5]) and the fact that 0 < g 
The denominator oj can be estimated as follows: 


tj-i < l/N. 


g / , \2 ^J—1 


N 

C\/e 


N 


> 


+ {a - tj-i) {tj - a) 

C\fs 

since (g - tj_i) [tj - g) > 0. 


Therefore, 

U-hj ^ Ce N 

hjH ~ aqN C-y/e aq 

This completes the proof of (HD). □ 

It is easy to see that Ajv is an L-matrix. The next lemma shows that is 
an M-matrix and that the modified discretization m is stable uniformly in e. 


Lemma 2. Let e be sufficiently small, independently of N, and let a > 4//3. 
Then the matrix An of the system (1101) satisfies 



< C. 


Proof. We want to construct a vector v = [r^o, vi,..., vn]'^ such that 

(a) Vi > S, i = 0,1,..., N, where <5 is a positive constant independent of both 
£ and N, 

(b) Vi < C, i = 0,1,... ,N, 

(c) cTi := kvi-i + diVi + ViVi+i > 6, i = 1,2,... ,N - 1. 

Then, according to the M-criterion, 

Ilia'll < <5-1 |H|<C. 
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The following choice of the vector v is motivated by HHUi: 


Vi = 


a — Hi + A, 
a — Hi + 


1 + PJ 


(1 + P) 


J—i 


i<J-l, 
i > J, 


where pj = [ihj/{2e), p = j3H/{2e), and a and A are hxed positive constants. 
Since HN < C, there exists a constant a such that Vi > a — Hi > S > 0, so the 
condition (a) is satisfied. Then, because of Vi < a + A, the condition (b) holds 
true if we show that X < C. We do this next as we verify the condition (c). 
When 1 < i < J — 2, we use (ED to get 


cTi = {k + dj + n)vi + kH - nH 


hi 

= JjCiVi - 

s 

E ^ hihi 

hi 

hi+1 hi+i 

IV 

1 

1 

e 

\ bi bihi 

hi+i 

J 2 2h,+i 


- hi) 

bi bihi 

hihi^i 

2 2hi+i 

>-l + h 

- a 2 

- 2 

- - d > 0. 
a 


The constant <5 exists because of the assumption a > 4//3. 
For j = J — 1, we have 


CTj-i 


hj. 


-cj-ivj-i + Ij-iH - rj^iH 


H 

+ ^ij-i + Xdj-i + r j-i 
e 


A 


> -- 


> -- 


hj-i 

e 

ij-i 


1 + PJ 

£ bj-ihj-i 


hj 

bj-i Xpj 

-’’J-i 


- rj-i 


1 + PJ 


Apj 

1 + PJ 


> g +-+ ( ^ 

- hj.i 2 \hjH hjH J2s + I3hj 
_ e B /2£ + + hj)\ XjBhj 


hj—i 2 Y 2hjH 



^j-i 


AH - 2 


J 2e -\- jBhj 


with a suitable positive constant A. We can choose such A because the estimates 
H < 2N~^ and q — tj_i < q — tj -2 < 1 imply 


Iff 


A/3 N 


TH - I-= 7T7 - — (9 - ^J-i) (9 - tj-2) >N[—-—]>0. 


hj- 


4ff 


aq 


A/3 


1 


aq 


























For i = J, we get 


aj = ^cjvj + IjH -rjH + X 

li 

> __L + — + 

- hj^ H 


'h + ^ + 


rj 


1 + pj (l + pj)(l + p)_ 


A 


(l + pj)(l + p) 

- H hj 2 

X 

^ (1 +Pj)(l + P) 
e{H-hj) 

- 2 hjH 


[^j(l + Pj)( 1 + p) + rfj(l + p) + r j] 


[^j(l + Pj)(1 + p) + rfj(l + p) + rj] 


> 5 > 0. 


The above estimate holds true because (IT^ implies that 

e{H-hj) ^ ^ 2 

hjH ~ aq ~ a’ 

when e is sufficiently small, and because we can show that 
[^j(l + Pj)(1 + p) + dj{l + p) + rj] > 0. 

Indeed, 

^,/(l + Pj)(l + p) + dj{l + p) + rj = Ijpj + Ijpjp - Tjp 

£ phj £ fihj pH 
^ ~ hjH^ ~ hjH^^ 
PH 


+ 


e bjhj 

P^ , Pbjhj 
4e 


2e 


2H£ 

P^ Phj{h.j + H) 
4£ 4H£ 

4e 4e 


Finally, when J+l<i<7V — 1, we have 


o-j = + kH - r^H - 


h 


1 + PJ 


A 


+ 


A 


(1 + p)® 1 

A 


(1 + p) 


i — J 


1+ PJ [(1+ p)*+i-‘^ (i + p)*-"' 




> - 
- 2 


(l+pj)(l + p)*+i-‘^ 

{k -ri+ kp)p 
{i + pj)ii + py+^-j' 
h p \ xpil+py-^-^ 


A 


H 2H 


1 + PJ 
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□ 

By examining the elements of the matrix Ajv, we see that 

\\An\\ < CN\ 

When we combined this with Lemma [U we get the following result. 

Theorem 3. The matrix An of the system (1101) satisfies 

k{An) < CN^. 


6 Uniform convergence 


Let Tj, i = 1,2,..., N — 1, be the consistency error of the finite-difference oper¬ 
ator , 

Ti = - fi. 


We have 


n = Ti[u] := - (Cu), 

and by Taylor’s expansion we get that 

\n[u]\<Ch,+,ie\\u'''\U + \\u"\U), (13) 

where H^Hi := l<?(3^)| for any C'(/)-function g. Let us define 


J ^TiW], ^ <i< J, 

[ Ti[u\, J+l<i<N — 1. 


Lemma 3. The following estimate holds true for all i = 1,2,N — 1: 

\fi[u]\<CN-\ 

Proof. We use the decomposition @ and estimates For the smooth part of 
the solution, it is easy to show that |f [s]| < CN~^. Then we need to show that 

\h[y]\<CN-\ 


Case 1. Let i > J + 1, i.e. ti_i >tj>a. Then we have 
\n[y]\ = \TM<Ch.+i (£|l2/"'ll. + ||y"ll.) 

< CN-^X'{U+i)e-^e-I^AU-i)/e 

< CN-^, 


where we have used the fact that e ^e °-h°‘ACV^ < (j 
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Case 2. Let i < J, i.e. ti_i < a, and at the same time, let ti-i < q — 3/N. 
Note that, when ti-i < q — 5/N, we have 

ti+i<q-l/N<q and q - ti+i > ^{q - ti^i). 

This is because 

3 2 2 

9 ~ ti-l — ^ ^ ~ ^i-l) — 

which gives 

2 1 2 2 1 

g - ti+i =q- ti_i - — = -{q- U-i) + -{q - U-i) - — > -{q - 

Therefore, 

< ^CK+i {e\\y"'\\i + \\y''\U) 

< CN-^ 

< CN-^ 

< Ce-^N-^{q - 

< CN-\q - 

< CN-^, 


because (g — < C. 

Case 3. In the last case, we consider the remaining possibility, g — 3/iV < 
ti-i < a. We use the fact that Cy = 0 to work with 


= -^Inl 


< [Pi + Qi + Ri) ■ 


where 

Pi — s\D yi\^ Qi — bi\D yi\^ and Ri — Cj||. 

We now follow closely the technique in m Lemma 5], (see also [UllH]), to get 


-[7 [Pi + Qi + Ri) ^ C 

£1 


5 (y- 


< CN-^. 


n 



_|_ g-/3'>'(‘i)/e 


□ 

Remark 1. The technique used in the above proof is based on [S], where 
the same approach is successfully applied to reaction-diffusion problems. This 
approach is originally due to Bakhvalov [I]. The technique works here for 
convection-diffusion problems © because an extra e-factor is obtained from 
the preconditioner (jH]). 
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When Lemmas [5] and [3] are combined, which amounts to the use of the 
consistency-stability principle, we obtain the following result. 

Theorem 4. Let s be sufficiently small, independently of N, and let a > 4//3. 
Then the solution of the discrete problem © on the VB-mesh satisfies 

< CN-\ 

where u is the solution of the continuous problem ©• 
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